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ABSTRACT 

In this study, a model for the tuberculosis infection 
considering vaccination and lost-sight compartement is 
formulated, there are six populations in this model, 
Susceptibled, vaccinated, exposed, lost sight, infected, and 
recovered. The lost sight populations are infected but do not 
get any treatment and still can spread the tuberculosis, the 
infected population are infected but already got a treatment 
and no longer spread the tuberculosis. The local stability are 
obtained by analyzing the epidemic threshold ,7? 0 . The result 
shows that the disease-free equilibrium is locally 
asymptotically stable when the condition 0l o < 1 is satisfied, and 
the unique endemic equilibrium exist and it is locally 
asymptotically stable if K 0 >1 is satisfied. The numerical 
simulation are also performed to support the analytical result. 
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I. INTRODUCTION 

Tuberculosis (TB) is one of the most infectious 
diseases caused by mycobacterium tuberculosis. It can 
attack various organs, mainly the lungs. This disease can 
lead to complications and death if the sufferer does not take 
any treatment or involved in incomplete treatment. 
Tuberculosis had allegedly appeared since 5000 years BC 
but the cure was just discovered in the last two century [7]. 

TB bacteria are spread through the air when 
infected person coughs or speaks. When a person nearby 
breaths the air, the bacteria will settle in the lungs and they 
can move to other organs. Individuals with TB will 
experience several symptoms such as bad coughs, weight 
loss and night sweating. The symptoms usually appear 
slowly within a few months and it impacts on the delay of 
the treatment. Without intensive treatment, two out of three 
cases of TB disease will lead to death [11], 

Worldwide, TB is one of the major public health 
issues. It is among one of the top 10 causes of death and the 
leading cause from a single infectious agent (above 
HIV/AIDS). In 2017 TB caused an estimated 1,3 million 
deaths among HIV negative people and also 300.000 deaths 
from TB among HIV-positive people. Globally, there are 10 
million people developed TB disease in 2017: 5,8 million 


men, 3,2 million women, and 1 million children [11]. From 
the latest data, Indonesia ranks the third place after India 
and China in terms of the infected number of TB in the 
world. There are estimated more than 1 million cases of TB 
in Indonesia but only 420.000 cases were reported to 
Ministry of Health [12]. 

The controlling of the TB cases in Indonesia were 
done by several strategies. Started with Directly Observed 
Treatment Short-Course (DOTS) at 1995-2005, followed by 
stop TB strategy at 2006-2015. Then according to the 
results of the prevalence survey, the strategy was revised to 
the TB elimination strategy. According to this strategy, the 
notification of tuberculosis cases increase drastically and 
the success of national tuberculosis treatment remains high 
at 87% [12], 

One of the challenge of the TB prevention 
program is Multi Drug Resistance (MDR) TB. This is the 
state when the body is resistant to at least izonaid and 
rifampin, the two most potent TB drugs [3], It worsens the 
patient’s condition and increases the infectious rate. It is 
caused by misuse or mismanagement of TB drugs treatment 
such as the wrong dose of drugs or incomplete treatment. It 
is more common in people who do not take their drugs 
regularly and continuously by reason of the cost of the 
drugs, the access to public health facilities, and side effect 
of the drugs. Reported cases of MDR-TB in Indonesia keep 
increasing since 2011 and not all of them were confirmed to 
retake the treatment [12]. 

In eradicating the infectious diseases including 
tuberculosis, one essential thing that should be considered 
is controlling the number of infection. Mathematics had 
been an important tool in analyzing the transmission 
behavior of the disease [5]. Many studies about 
mathematical models and its analysis of TB infection had 
been conducted. Temgouaef al. [8] developed the model for 
tuberculosis with lost sight and multi-latent compartment. 
Another study conducted by Aprilianief al. [1] analyzed the 
SEIR model along with the vaccination compartment. This 
study combines the existing models mentioned before to 
analyze the effect of treatment, treatment adherence, and 
vaccination in suppressing the spread of TB infectious. 

II. MATHEMATICAL MODEL 
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The total population denoted by (V(t) is divided 
into seven compartments, S(t) represents the susceptible 
population, V (t) represents the number of infected 
vaccinated population, E (t) represents the number of 
exposed population, L(t) represents the number of lost 
sight population, /(f) represents the number of infected 
population, and R(t) represents the number of recovered 
population. 

In this model, the susceptible individual will go to 
the exposed population first before becoming infectious. 
The exposed individual will progress either to infected 
population or lost sight population. Individuals at Infected 
population are given a treatment and progress to recovered 
population after they complete the treatment. Lost sight 
population are the individuals from exposed population 
who does not take any treatment and individuals from 
infected population who does not complete the treatment. It 
can also go back to the infected population whenever they 
retake the treatment. We consider that the infected 
individuals in treatment cannot spread the disease because 
they are under surveillance, therefore the disease can only 
be spread by lost sight population. Vaccinated populations 
can also be infected but in lower rate depending on the 
vaccination effectiveness rate. The flow diagram is shown 
in Figure 1. 



Fig 1. Diagram flow of tuberculosis infectious 


Based on the diagram flow above, we formulate 
the following system of ordinary differential equation: 
dS 

— = V - fSL - ( M + B)S, 

J = 6S-iiV-{l-a)pVL, (1) 

dE 

— = (1 - p)PSL + (1 - a)pVL - (jj. + v + c)E, 
dt 

dL 

— = pfESL + qvE + 0(1 - c)/ - (p + + c)L - 5(1 

dt 

— c)L, 


— = p(l - f)BSL + (1 - q)vE + 5(1 - c)L - (jj. + p t + c) 
dt 

-0(1 - c)/, 
dR 

—— — c(E + /, + /) — pR, 
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by: 


III. RESULT AND DISCUSSION 

System (1) has a disease-free equilibrium given 


T 0 (S,V,E,L,I,R ) = (— ,■ 


Y e 


0 , 0 , 0 , 0 ), 


B+p n{6+ii) 

And an endemic equilibrium given by: 

T*{S,V ,E,L,l,R) = (, S*,V*,E*,L*,r,R *), 

Where 
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Basic reproduction number, denoted by JIq, is an 
expected number of secondary cases produced by a typical 
infective individual in a completely susceptible population. 
It is obtained by the next generation matrix formulated in 
[4]. We obtain the basic reproduction number as follows: 


72, 


where 
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Theorem 1: The disease-free equilibrium is locally 
asymptotically stable if Jig < 1. 

Proof: After linearization for system (1) at Tg we 
have a Jacobian matrix as follows: 
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where: 
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The eigen values are obtained by |/ ro — 2i| = 0 
therefore we have the characteristic equation as follows: 

(4 - 7n)(4 -Jn)^ — 766)C-4 3 + a i4 z + a 2 A + a 3 ) 

= 0 

where 
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From equation (3) we obtain the eigenvalues as 

follows: 

4-1 = -0 - p, A 3 = -p A 2 = -p 

All the eigenvalues above are negative, the others 
will be obtained by solving the following equation: 

A 3 T ct x A 2 + cl 2 A T cl 3 (4) 

There are three conditions that must be satisfied 
for the cubic equation above: 

A4 + 2.5 + A (, — — cl x 
A 4 A s + A 4 A 6 + A s A 6 = a 2 (5) 

2 4 2 5 2 6 = —a 3 

Let 72 0 < 1, such that according to equation (2) we 

have: 

72g + 72(j + 72$ + 72$ 1 ( 6 ) 

Since 72$ > 0, 72$ > 0, 72$ > 0, 72$ > 0,therefore 

72$ < 1, 72$ < 1, 72$ <1 , 72$ < 1. 72 1 < 1 gives: 
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Then 


Inequalities ( 8 ) and (9) gives a 2 > 0 and we also 
assume that 72$ (—-—-) < lwhich gives a 3 > 0 . 

V 0 ( 1 — G)J 

Therefore, the conditions (7) will be satisfied when 2 5 < 
Oand A 6 < 0. 

Since all the eigenvalues are negative when 
72 0 < 1, according to [9], it is proved that the disease-free 
equilibrium is locally asymptotically stable if 72 0 < 1 . 


Theorem 3: the endemic equilibrium is locally 
asymptotically stable if 72 0 > 1 


Proof: According to Castillo-Chaves and Song [2], let 
(p = P 3 be a bifurcation parameter. When 72 0 = 1, we have 

<P = V‘ 

n(6 + /i)(c+/i + (it)(c +ii + v)(J + n + /it + ip - c(—1 + 5 + 0)) 
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The disease-free equilibrium has five negative 
eigenvalue and one zero eigenvalue. The zero eigenvalue 
has a right-side eigenvector (u 1 ,u 2 ,u 3 ,u 4 ,u 5 ,u 6 ) and a 
left-side eigenvector (v 4 , v 2 , v 3 , v 4 , v 5 , v 6 ) where 
u x < 0 , u 2 < 0 , u 3 > 0 , u 4 > 0 , u 5 > 0 , u 6 > 0 , v 4 — 0 , 
v 2 = 0 , v 3 > 0 , v 4 > 0 , v 5 > 0 , andv 6 = 0 


Now we define 
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d 2 f 2 

dx x dx 4 

d 2 fs 


(T 0 ,<pl = (1 - p)<P* 
(T 0 ,(pl = p(i -/>* 


dx 1 dx 4 
also partial derivatives for b: 

q 2 -C 
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IV. NUMERICAL RESULT 

We perform the numerical simulation to support 
the analytical result. The simulation visualizes the 
dynamical population in both condition, Jl 0 < 1 and 


d 2 fs 

dx 4 dcp 

d 2 f 4 


(T 0 ,(pl = P(1 ~f)x 1 


C T 0 ,(pl = Pfx 4 


dx 4 d(p 

According to equation (12) we obtain a < Oand 


b > 0. 


Consequently, when (p changes from (p < cp *to 
(p > (p*, the free-disease equilibrium changes from stable 
to unstable and the endemic equilibrium changes from 
negative to 

positive and it is locally asymptotically stable. It proves that 
the endemic equilibrium is locally asymptotically stable if 
72 0 > 1 . 


72 0 > 1 according to theorem 1 and theorem 2. The initial 
conditions we use for this simulation are: 5(0) = 
10.000, T(0) = 50.000, £(0) = 10.000,L(0) = 5000dan 
P(0) = 0. And the parameter values are presented as 
follows: 


Table 1 Parameter Values 


Parameter 

Value 


Source 

72 n < 1 

> 1 

V 

1 

1 

[1] 

P 

-0.996691 

-1.06243 

[1] 

Pr 

-0.302297 

-0.148983 

[1] 

P 

0.236319 

0.578424 

[1] 

9 

1 

1 

[1] 

o 

0.10979 

0.127638 

[1] 

c 

0.139156 

0.129349 

[1] 

V 

-0.89669 

-0.052261 

[1] 

s 

0.375597 

0.0997238 

[8] 


0.73202 

0.389793 

[8] 

p 

-0.603326 

-0.648159 

[1] 

f 

-0.778148 

-0.553713 

[1] 

q 

-0.187724 

-0.107216 

[1] 
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Figure 2: the population dynamic at 72 0 < 1 an d > 1 


V. CONCLUSION 

We’ve formulated a model for TB transmission by 
combining the existing model with vaccinated and lost sight 
in a form of ordinary differential equation system. The 
system gives two equilibrium, a free-disease equilibrium 
and an endemic equilibrium. A free-disease equilibrium is a 
state when there is no outbreak in population, and endemic 
equilibrium is a state when the disease is endemic in 


population. 

The stability analysis shows that the disease-free 
equilibrium is locally asymptotically stable if < 1 is 
satisfied and the endemic equilibrium is locally 
asymptotically stable if > 1 is satisfied. The numerical 
simulation shows the suitable result with the stability 
analysis. 
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